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Abstract 

The use of linear response theory for forced dissipative stochastic dynamical systems 
through the fluctuation dissipation theorem is an attractive way to study climate change 
systematically among other applications. Here, a mathematically rigorous justification 
of linear response theory for forced dissipative stochastic dynamical systems is de- 
veloped. The main results are formulated in an abstract setting and apply to suitable 
systems, in finite and infinite dimensions, that are of interest in climate change science 
and other applications. 

1 Introduction 

One of the cornerstones of modern statistical physics is the fluctuation dissipation the- 
orem (FDT) which roughly states that for systems of identical particles in statistical 
equilibrium, the average response to small external perturbations can be calculated 
through the knowledge of suitable correlation functions of the unperturbed statistical 
system. For some of the many practical applications of the FDT, see for example 
IIKTH85IIBal97l . The low frequency response to changes in external forcing for var- 
ious components of the climate system is a central problem of contemporary climate 
change science. Leith IILei75l suggested that if the climate system satisfied a suit- 
able FDT, then the climate response to small external forcing should be calculated by 
estimating suitable statistics in the present climate. The climate system is a forced 
dissipative chaotic dynamical system which is physically quite far from the classical 
physicists' setting for FDT. Leith's suggestion has created a lot of recent activity in 
generating new theoretical formulations and approximate algorithms for FDT that were 
apphed to chmate response with considerable skill ( IIGD99I IGBD021 IGB07I IGBM08I 
IAM07I lAMOSl I AM09I IMWTOl ) . 

The goal here is to provide a rigorous justification of linear response theory at in- 
finite times for a rather large class of general stochastic systems, including a class of 
problems with potential applications to climate change science. The improved regu- 
larity due to the presence of noise simplifies key aspects of the problem considerably 
and circumvents the difficulties found to occur even in relatively simple deterministic 
systems (see ||Rue97| iBaMl for an overview and HE Y93I |KL99l IGL06I iGLOSl iBSOSl 
for a number of advances in this direction). As a consequence, we are able to focus on 
systems that have some other features of realistic systems, namely a lack of ellipticity 
(we will only require hypoellipticity), non-compactness of state space, a lack of global 
Lipschitz continuity of the coefficients, and even infinite-dimensionality of the state 
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space. The results established below apply to a very broad class of nonlinear function- 
als of state space which trivially include many quantities of interest, such as the mean 
and variance of subsets of variables. 

The mathematical formulation of FDT as linear response theory for forced dis- 
sipative stochastic dynamical systems IIRis89l IMAG05I IMFK08I IPalOll IMWIOI is 
an appropriate setting for these applications since many current improvements in the 
comprehensive computer models for climate change involve stochastic components 
IPWIOI , while lower dimensional reduced models typically also involve stochastic 
noise terms I MFK08 . IM AG 1 01 IMG Y 1 01 . While there is a recent formal systematic 
discussion of FDT for time-dependent stochastic systems of interest in climate change 
science IMWIOI . it is not obvious whether the results of this article generalise in a 
straightforward way. The remainder of this article is structured as follows. After in- 
troducing our setup and notations, a general abstract formulation of the main Theorem 
is presented in Section |2] with a general discussion on its applicability provided in 
Section [3] Section |4] illustrates the scope of the theorem for finite-dimensional forced 
dissipative stochastic systems with a prototype structure which is relevant for some 
applications in climate change science. 
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2 Main result 

Throughout the entire article, we will consider a family of Markov evolutions {Pf : 
t > 0, a G R} on a separable (either finite or infinite-dimensional) Banach space B. In 
other words, the are Markov operators over B such that Pq(p = tp and VfT^f = 
Vg^^ip for every bounded measurable function — > R. In this abstract setting, we 
take the pragmatic view that a Markov operator V is determined by a transition kernel 
Pas 



where the map A ^ P{x, A) yields a Borel probability measure for every x G B and 
that the map x P{x, A) is measurable for every Borel set A. We will denote by V* 
the adjoint of V acting on the space of finite signed Borel measures on B by 



We will also fix a particular value oq G R and we are interested in the long-time 
behaviour of the system described by "P" for a close to oq. Throughout this article, we 
assume that for every a, there exists a probability measure yu^ on B which is invariant 
for Vf, so that {VI;)* iia = /^a for every t > 0. The aim of this note is to show 
that, under very weak natural assumptions on the dynamic, the map a n- {ip,ij,a) is 
differentiable at qq for every sufficiently regular test function (p: B R. (Note that 
the fact that a G R is not a restriction since if we take a G R" instead, it suffices to 
consider its components one by one, keeping the others fixed.) 
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One of the main features of our proof will be that it relies on rather 'soft' arguments 
and has conditions that are easily verifiable in practice. For example, it does not require 
the transition probabilities to have a density with respect to some reference measure 
(typically Lebesgue in finite dimensions or Gaussian in infinite dimensions). 

Denote by C^(B) the set of all functions ip: B ^ R such that there exists > 0, 
a linear map T: B ^ R^, and a smooth, compactly supported map ip : — > R such 
that (fi = (f o T. Given two continuous functions V,W : B [1, oo), we set Cy to 
be the closure of C^{B) under the norm 

11.11... = ™p(^ + ^). (2.1) 

x^B^ V{x) W{x) J 

If we quotient this space by the space of all constant functions, then this norm is equiv- 
alent to the Lipschitz norm corresponding to the distance function 

dvMx, y) = Qw{x, y) A iV(x) + V(y)) . (2.2) 

Here, for any positive continuous function W bounded away from 0, we define the 
metric ^ 

gw{x,y)= inf / Wi-/(s))\\^(s)\\ds , 

where the infimum runs over all differentiable curves 7 : [0,1] B with the prescribed 
boundary conditions. (Note that if = 1, then one simply has Qwix, y) = \\x — y\\.) 
We then make the following assumption: 

Assumption 1 There exists a time t > and a constant g < 1 such that ||'P°°(/? — 
{•^,i^ao)\\i;V.W < q\W- ^^aa)\\l■y,w for every ip S Cy^w 

Remark 2.1 This spectral gap condition is quite different from the or L'p spectral 
gap conditions often encountered in the literature. It is however possible to verify that 
such a condition is satisfied for a very large class of finite-dimensional and infinite- 
dimensional (hypo-)elliptic diffusions. See II HMOS I for the verification of Assump- 
tion[T]in the case of the 2D stochastic Navier-Stokes equations and Section|4]below for 
a class of finite-dimensional SDEs. 

In particular. Assumption [T] implies that the invariant measure for V^" is unique, 
that the spectrum of as an operator on Cy has an isolated eigenvalue of multi- 
plicity 1 at 1, and that the remainder of its spectrum is contained in the ball of radius 
Q around the origin. Therefore, for every centred function G Cy there exists a 
unique centred function -0 G Cy [y such that 1/' — Vi°i\j = ip. We will henceforth use 
the notation (/' = (!- 7'"")" V- 

Our second assumption concerns the regularity of Vi with respect to the parameter 
a. In order to state this assumption, we introduce the space Cjj, which is the weighted 
space of continuous functions obtained by completing C^{B) under the norm 

II II \f(x)\ 
Mu = sup—— . 

x^B U{X) 



We then assume that 
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Assumption 2 There exists a continuous function U > V from B to [\, oo) such that, 
for some fixed i > and every ip S Cyy/, the map a H- Vf^ is differentiable in a 
neighbourhood of gq when viewed as a map from R to Cu- Denoting this derivative by 
P'^ip, we furthermore assume that 



for some constant C > 0. 

Remark 2.2 Note how we are allowed to lose one derivative as well as some growth 
in the norm of ip when computing P^^s. This is what makes the verification of ( 12.3b 
very easy in practice, provided that U, V, and W are chosen in a judicious way, see 
Section Hlbelow. When considering such problems on spaces of probability measures 
instead of spaces of test functions, it has been know for some time that in the context 
of dynamical systems, spaces of functions were suitable to take advantage of the 
expanding character of a map ||Rue89| . Several highly non-trivial extensions of this 
fact were considered much more recently in IIGL06I IbTOTI IBS08 I . See also IIAGY06I 
for an instance of the use of weights in a related dynamical systems context. 

Finally, we assume that we have an a priori bound on the integrability of the in- 
variant measures: 

Assumption 3 There exists e > such that supj^_^^j^^ U (x) fia(dx) < oo. 
Our main result then states that: 

Theorem 2.3 Let {7'f}aeR be a family of Markov semigroups over a separable Ba- 
nach space B such that there exist functions U,V,W : B R+ such that Assump- 
tions\l]^and\3}are satisfied for some a — and some fixed t > 0. Then the map 
a I— {(p, ^a) '■s differentiable at oq for every tp £ Cy ^ and the identity 



holds. In particular, the right hand side of H2.4\l is well-defined. 

Remark 2.4 As pointed out by one of the referees, the spectral gap condition is not 
essential for Theorem l2.3l to hold. It could be replaced by the requirement that the map 
a I— >■ /Xo is continuous in the topology of weak convergence and that the test function 
ip is such that there exists a function ij.' G Cy^y such that (1 — 'Pt")i' ~ V — {f, l^aa)- 
It is not clear how easy the latter condition would prove to verify in practice for 
systems that do not possess a spectral gap as it would require a detailed understanding 
of the properties of the resolvent which is bypassed by Assumption[T] 

Remark 2.5 The formula (12.4b is identical to that familiar to dynamicists, see for ex- 
ample BBalOSI Eq. 3]. Since the left hand side in (12.41 ) is independent of t, we can take 
the limit t — > in the right hand side. If the test function (p is centred with respect to 
fiao ™d sufficiently 'nice' so that t(l — Vf")^^(p ^aof' we then obtain the linear 
response formula 



Ilf; 



(2.3) 




(2.4) 




(Lao ^aof^l^ao) , 



(2.5) 
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where is the formal derivative of the generator La with respect to the parameter 
a. (This can also be obtained by formally differentiating the identity >C*/ia = with 
respect to the parameter a.) 

If the state space is R", jUq^ has a density exp(— 7?(a;)) with respect to Lebesgue 
measure and L^^ = F(x)'Vx is a first-order differential operator, then we can define the 
conjugate current J by 

J{x) = F{x)VxH{x) - div, F{x) . 
With this notation, ( 12.5b can be interpreted as the Green-Kubo relations 



-^{ip^Ha) = {JC^^ip, flao) ^ I ^{J{XQ)ipiXt))dt , 

dCL a—ao 



where expectations are taken over the stationary Markov process xt with invariant mea- 
sure flao- 

Remark 2.6 Note that even though we assume uniqueness of the invariant measure for 
Vf" , this need not be the case for a uq. Even though uniqueness doesn' t need to hold 
for a ao, our result implies differentiability of the map a ((^, /ia) at a = ao for 
any choice of a M- fj,a- This is analogous to the result shown for partially hyperbolic 
systems in IIDol04ll . 

Remark 2.7 The classical perturbation results found for example in IIKat76l do not 
seem to be applicable directly to our situation since we do not assume that the map 
a ^ P^ip is differentiable in Cy 

Proof. Let us first show that the map a Ha is Lipschitz continuous at a = ao if 
we endow the space of probability measures on B which integrate dy.wiO, x) with the 
norm dual to that of Cy^^. For every if G Cy ^, we have the identity 

{ip, Ma - Mao) = {Vt\, Ma - Mao) + {(V^ - V^V', Ma) . 

Since none of these expressions changes if we add a constant to ip, we can assume 
in the sequel that ip is centred with respect to i^iag so that, using the invertibility of 
(1 - Vt") and setting ip ^ (1 - Pt°)^- we obtain the identity 

ii;, ila - ilao) - {{Vt - ^r)(l - n"")" V, Ma) - (2.6) 

It now follows readily from Assumptions |2] and |3]that 

sup {lp,^^a-^ia,)= sup ((P,'' - 7^ )(1 ^ ^D' V, Ma) 

||'/'l|l;V,W = l ||'P||l;V',W = l 

<C\a-ao\ sup Ma) < Cja - ao| , 

H'\\u=i 

for some positive constants C and C, as required. 

Since the norm dual to that of Cy ^y (restricted to differences between probability 
measures) is nothing but the Wasserstein-1 distance corresponding to dv,w IIVil03l . it 
follows that — > Mao weakly for a oq. Combining this with Assumption |3] we 
conclude in particular that 

(*,Ma) ^ (*,Mao), V*GCc/. (2.7) 
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(This follows by a standard 'diagonal argument' from the fact that '5 can be approxi- 
mated by continuous bounded functions by the definition ofCjj-) 

We now have all the preliminaries in place to show that {ip, jia) is differentiable 
and to compute its derivative. It follows from Assumption |2] that, for every (centred) 
S Cy^, there exists a function a i~> TZa G Cij with lima^Q,-, ||7?.a||(7 = and such 
that we have the identity 

^ -{vt - v^"){i - v^°r^^ = p°°(i - T'ro^V + T^a ■ 



a — flQ 

As a consequence of ( 12.6b . we have 

+ {(P^ ~ rr'){l ~ TT')- V, Ma - Mao) , 

SO that 

^^'^"'^""^ - (p^d -^r)-V,Mao) < 2|(7^,,A.,„)| + |(7^„,A^„)| 

a — ao 

+ |(P,"«(l-7'r)~V,Ma-Mao)l ■ 

The first two terms on the right hand side converge to zero thanks to the continuity of 
TZa at ao combined with Assumption[3j while the last term converges to zero thanks to 
( 12.7b . thus concluding the proof. □ 



3 Domain of applicability 

In this section, we discuss under which circumstances one can expect a spectral gap 
result in a norm of the type ( 12.1b . Harris' theorem IIHM09bllMT93l provides a general 
methodology to verify whether a system satisfies a spectral gap in Cy for some V . If it 
so happens that such a system also has the property that one can find W such that, for 
some fixed to > 0, is a bounded operator from Cy into Cyy^, then it follows that 
Vl° has a spectral gap in Cy ^. However, such an approach presents two problems: 

1. While Harris' theorem is suitable for finite-dimensional diffusions, there are 
many examples of stochastic PDEs satisfying a spectral gap in a space of the 
type Cy for which Harris' theorem does not apply. 

2. In finite dimensions, as soon as a system satisfies Hormander's bracket condition, 

maps continuous functions into smooth functions. In particular, there always 
exists some W such that it maps Cy into Cyy^- The problem is that there does 
not appear to be a general result available that controls the growth of such a 
function W . This is problematic, since the function I] appearing in Assumption|2] 
is generally comparable to W , so that Assumption [3] is hard to verify without 
control on W . 

Our aim is therefore to present conditions under which a spectral gap in a space of 
the type Cy ^y can be verified directly, without requiring a spectral gap result in a space 
of the type Cy. In IIHMS09I . a general framework was laid out that gives conditions 
under which a "spectral gap" is obtained in a Wasserstein-type distance. Unfortunately, 
the Wasserstein distance in question is with respect to a metric comparable to the square 
root of the original metric on the space. Therefore, in our setting, this actually yields 
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a spectral gap in a space of i-Holder continuous functions, which is useless for our 
purpose. On the other hand, a spectral gap in a space of the type ( 12.1b was obtained for 
the stochastic Navier-Stokes equations in |HM08J. 

In this section, we obtain a slightly different set of conditions under which a spectral 
gap in a norm of the type ( 12.1b can be shown for a general Markov operator V over a 
Banach space B. Our conditions are inspired by those in [HMOS, HMS09 |, aiming in 
particular at being applicable to equations with conservative quadratic nonlinearities, as 
arising in applications. The main ingredient is again a gradient bound of the following 
type, which was shown to hold for a large class of stochastic PDEs in ||;HM06. HTs409al : 

Assumption 4 There exist continuous functions Ui,U2'- B ^ such that, for every 
e > there is a constant > such that the bound 

\\Drip{x)\\^ < eU!(x){r\\D^\\^)(x) + C,Ui(x){Pip^)(x) , (3.1) 

holds for every x ^ B and every Lp G C^iB). 

Remark 3.1 In the case of the stochastic 2D Navier-Stokes equations, this bound can 
be shown to hold with Ui{x) = U2{x) = exp(77||a:|p) for arbitrarily small 77 > 0, with 
x the vorticity of the vector field and its L^-norm. 

The aim of this section is to show that if a Markov operator P satisfies Assump- 
tion |4] and has sufficiently good contraction properties, then it is possible to deduce 
a spectral gap in the norm ( 12.11 ) (for a suitable choice of V and W), provided that 
the transition probabilities satisfy a kind of topological irreducibility condition. This 
statement can be formulated precisely in the following way: 

Theorem 3.2 Let V be a Markov operator over a separable Banach space B mapping 
Cq^(B) into Cy-:^ and satisfying Assumption^ Suppose furthermore that there exist 
continuous functions V,W : B R-|- satisfying 

U^PW^ + UiPV^ < CW^ , PV <^V + K , (3.2) 

for some constants C and K. 

Finally, assume that there exists a point x^, € B such that, for every e > and 
every C > there exists a > such that 

inf P(x,B,(x^))> a . (3.3) 

X : V(x)<C 

Then, P has exactly one invariant probability measure fii,. Furthermore, there exist 
constants C and 7 > such that the bound 

lirV - ^i.)\\l■,v,w < Ce-^"\\^ - fi.)\\i,v,w , (3.4) 

holds for every ip G Cy yy. Finally, there exist constants 5 > 0, (3 > and g < 1 such 
that the bound 

\\Pip - (<y5, M*>l|l:l+/3V,6-iW < Q\W " ((^7 A^*) I! l;l+/3V,6-i W , (3.5) 

holds for every ip ^Cyy^. 
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Remark 3.3 It will follow from the proof that it is sufficient that the second inequality 
in (13.2b holds with i replaced by any other constant k < 1. Furthermore, we do not 
need (13.3b to hold for every C > 0, but only for some C > K/(l — k). 

Proof. Since the norms || ■ \\i-^i+f3v,5-'^w ^"d || ■ are equivalent for every choice 

of 13,6 > 0, ( 13.4b follows immediately from (13.5b . so that it is sufficient to check (13.5b . 

The proof uses very similar ideas to the proofs of IIHMS09I Theorem 4.7] and 
0HMO8I Theorem 3.4]. In particular, we use the following trick. For ^ > and (3 > 0, 
we introduce the distance 

dpA^, y) = 5-^Qw{x, A (2 + l3V(x) + pV{y)) . 

This distance is of course equivalent to the distance dv,w introduced in (12.2b . but it 
turns out that allowing the freedom of choosing both S and /3 sufficiently small will 
considerably simplify the proofs. With this definition, a function if is Lipschitz 
continuous with Lipschitz constant 1 with respect to dp^s if and only if 

\\Dip{x)\\ < S~^W{x) , \ip{x)\ < 1 + pVix) . (3.6) 

(For the second inequality, one might have to add a suitable constant to ip.) Denote by 
Lip^- p the set of all such functions. As in IIHMS09I IHM08I . we show that it is possible 
to choose 6 and /? in such a way that the bound 

dpA^x, ■ ), Piy, ■ )) < adpAx, y) , (3.7) 

holds for some a < 1 uniformly over all pairs x,y ^ B. We again show (13.7b separately 
in three different cases and we use separately the three ingredients of the theorem in 
each of these cases. 

The case qw{x, y) < S{2 + f3V{x) + l3V{y)). In this case, we make use of 
the gradient bound ( 13.1b . together with the 'super-Lyapunov' structure (13.2b to deduce 
that if If satisfies (13.6b . then for every e > there exists such that the bound 

\\Drip(x)\\ < ed-^W{x) + CeW{x) , 

holds uniformly for all such ip and for all /3 < 1, say. It follows that by first choosing 
£ = \ and then choosing 5 small enough so that < 1/ (4(5), one has 

\\DVv{x)\\ < ^W(x), 
which immediately implies that 

dfs,s{V{x,-),P{y,-))< sup \V^ix) - Pipiy)\ 

< sup inf / WDVipijis))]] |7(s)| ds 

1 

< — supinf / W{j{s)) \j(s)\ ds 

< -^Qwix,y) < l-dp^six,y) , 



as requested. 
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The case Qw(x,y) > S{2 + f3V(x) + l3V(y)) andV{x) + V{y) > 4(K + 2). 
In this case, we only make use of the fact that is a Lyapunov function. We have 
indeed 

dpA'P(x,-),ny,-))< sup \r^(x)-Vv{y)\< sup {\Vv(x)\ + \r(p{y)\) 

< 2 + (3rV(x) + prViy) < 2 + ^{V(x) + V{y)) + I3K 

< (2 - 2/3) + ^{V{x) + V{y)) < (1 - (/? A \))dpAx, y) , 

which again yields a contraction, but with a strength that depends this time on the 
parameter 13. Finally, we have 

The case Qw(x,y) > S{2 + f3V(x) + f3V(y)) and V(x) + V(y) < 4(K + 2). 
In this case, we make use of our final assumption, namely ( 13.3b . At this stage we 
assume that (5 > is fixed, sufficiently small so that our first step goes through. We 
can then find some sufficiently small e > so that dp^Xi,^ y) < \ for all y G B^iXi,), 
uniformly over /3 < 1. In this case, we can decompose the Markov operator V into a 
combination V = aVi + (1 — q)V2 of Markov operators such that Vi{x, B^{XiS) = 1 
for every x such that V{x) < A(K + 2). We conclude that 

dp^Vix, ■ ), • )) < adp^Viix, ■ ), Vi{y, ■)) + {!- a)dp.s{V2{x, ■ ), ^2(2/, ■ )) 

< a + (1 - a){2 + PV2V{X) + /^Ps^W) 

< a + 2(1 - a) + pVVix) + pVV{x) 

< 2 - a + ^{V(x) + V(y) + 2K) <2-a + PCdK + 4) . 

We can now choose p sufficiently small so that this constant is strictly smaller than 2. 
Since on the other hand one has dp^x, y) > 2, the claim now follows. □ 

4 Application to stochastic differential equations 

In this section, we aim to apply the results of the previous sections to a concrete class 
of stochastic differential equations. In order to keep our conditions clean, we restrict 
ourselves to SDEs with additive noise and polynomial nonUnearities. We thus assume 
that for some integers N, M > 1, we have 

N M 

dx(t) = ^Afkix, ...,x)dt + ^ak dwk(t) , a;(0) G R'' . (4.1) 

fc=0 k=l 

Here, the functions Afk are symmetric multilinear maps of order k and the ak are ele- 
ments of R''. We also define the function A/": R'' R"^ by 

N 

AA(x) = ^A4(a;,...,x), 

as well as the d x M matrix E given by E = (di , . . . , (Jm)- 

When N = 2, equations of this type typically arise as effective dynamic for large- 
scale structures in climate models IIMAG05I . Fixing N and M, we make the following 
structural assumptions on the A4's and the o-fe's: 
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Assumption 5 There exist constants c,C > such that the bound 

(a;,A/'(a;)) <C-c||a;||", Vx G R'' (4.2) 

holds. 

Remark 4.1 The bound (14.2b can hold even if the degree N of the nonlinearity is 
even. In this case, it impHes that {x, J\fN{x, . . . ,x)) = 0. If is odd, then (14.21 ) would 
actually also hold with A^ replaced by A^ + 1, but this is not important for our purpose. 

Assumption 6 Define Ak C R*^ recursively by Aq = spanjcri, . . . , crj\/} and 

Ak+i = span(Afe U {A/'jvCyi, ■ ■ ■ ,2/7v) : Vj e ^fe}) ■ 
Then, one has Ad = R^. 

Assumption 7 The control system associated to ( 14.71 ) is approximately controllable. 
In other words, for every initial condition x^, every T > 0, every terminal condition 
XT, and every £ > 0, there exists a smooth control u £ C°°([0, T], R*^) such that the 
solution to 

x(t) = M(x{t)) + a u{t) , x{0) = xo , 
satisfies \\x(T) — xt\\ < £• 

Remark 4.2 If A^ is odd, then Assumption [T] is redundant since it then follows from 
Assumption|6] see IIJK85I . However, in the more interesting case where A^ is even, the 
situation is not so easy. A number of criteria are given in IISus87l |Kaw90l IMBB03L 
but they are all sufficient without being necessary. In the particular case of truncations 
of the two and three-dimensional Navier-Stokes equations where furthermore E is as- 
sumed to be diagonal in the Fourier basis, it was also shown that Assumption|2]follows 
from Assumption l6l EMO 1 1 lRorn04l I AS05 1 . 

For given M, N > 0, let ^ denote the set of (Af, S) such that Assumptions |5j 
|6] and Q hold. Note that (Af, S) can be viewed as a vector in R" with n = d{M + 
Ef=o('^fc''))' so that J is nothing but some possibly rather complicated open subset 
of R". With these assumptions in hand, we have the following result, which follows 
immediately from IIDPZ96||MT931 : 

Theorem 4.3 For every (Af, S) G the system 114. l\l possesses a unique invariant 
probability measure ijl(X,t,)- 

Proof. The existence of an invariant measure follows from Assumption [5] by Krylov- 
Bogoliubov. The uniqueness is then a consequence of the strong Feller property (which 
follows from Assumption |6]) and the topological irreducibility (which follows from 
Assumption |7]) of the associated Markov semigroup. □ 

The aim of this section is to prove the following result: 
Theorem 4.4 Fix M, N >! and let ip: K"^ Kbe a function such that 
sup {\(pix)\ + ||D(p(a;)||)e-*ll^ll' < oo , 

for every 6 > 0. Then, the map (Af, E) i~> J^^ 'P(x) p^{j\f .T.){dx) is on J" and its 
derivative is given by the formula ( 12.41 ). 
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Before we turn to the proof of Theorem 14.41 we start by providing the a priori 
bounds on the solutions to ( 14.1b required for our analysis. For this, we first introduce 
the Jacobian process of (14.1b . which is the d x d matrix-valued two-parameter process 
Js,t given by the solution to the following random ODE: 

^ = DM{x)Js.t , Js.s - Id . (4.3) 
dt 

We also introduce the second variation process jfl which is a bilinear map-valued 
process given by 

rfj^^Vf C) 

for every ^ and ^ in R''. We then have the following: 

Proposition 4.5 Let (Af, E) be such that Assurnption\5\holds. Then, there exist strictly 
positive constants rj, v and C and, for every e > and p > there exists a constant 
Cjj^s such that the bounds 



Eexp(77|ix(Qf + 77^ ||a;(s)||^ds) < Cexp(77e-''*||a;(0)|l2) 



(4.4) 



mJs,tr+\\jflr) < Cp,eexp(e||x(0)f), (4.5) 



hold for every initial condition x(0) G R"^, and for every s,t Cz [0, 1]. 

Proof. Applying Ito's formula to the function x i-t- ||a;|p, it follows from our assump- 
tion that 



1 

-d\\x{t)\\^ < -c\\x\\^ dt + C dt + J2{x{t),<Jk) dwk(t) 



k=l 

for a possibly different constant C. Note that this inequality should really be interpreted 
in integral form. The bound ( 14.41 ) then follows from the version of the exponential 
martingale inequality given in I1HM08I Lemma 5.1]. 

In order to obtain the bound ( 14.51 ). note that it follows from (14.3b that there exists a 
constant C such that 

\\Js,tr < exp(pC^*(l + llx(r)ir-i)dr) . 

The requested bound now follows immediately from ( 14.4b . noting that for every 6 > 
there exists a constant Cp,s such that the inequality 

pCil + \\x\\'^-')<5\\x\f + Cp,s, 

holds for every .t 6 R''. □ 

The main ingredient for all probabilistic proofs of regularising properties for SDEs 
is the Malliavin matrix of the process. This is the R'^^'^-valued matrix Mt defined by 

Mt= I Js,t^T.*J*ds , 



where Jg.t is the Jacobian of the flow as above and E is the noise matrix defined at the 
start of this section. Note that A^t is a random object which furthermore depends on 
the initial condition xq of ( 14.11 ). We have the following a priori bound: 
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Proposition 4.6 Let (A/", S) be such that Assumptions\5\and\6\hold. Then, for every 
t > 0, every 5 > Q, and every p > 1 there exists a constant C > such that the bound 

EllA^r^r <Cexp(5||xof), 

holds for every xq G R''. 

Proof. This is an immediate consequence of ||HM09a| Thm 6.7]. The assumptions of 
that theorem are satisfied by Proposition 14.51 if we choose ^l^oix) = exp(77||a;p) for a 
sufficiently small constant 77 > 0. □ 



As a consequence of Proposition |4j6] we obtain the following gradient bound for 
the Markov semigroup Vt generated by solutions to ( 14.1b : 

Proposition 4.7 Let (J\f, E) be such that Assumptions\5\and\6\hold. Then, for every 
t > and every rj > there exists constant C such that the bound 

\\DVMx)f < Cexpir,\\xW')iVt^^)(x), 

holds for every x G R'* and every Lp G C^(R'^). 

Proof. The proof is almost identical to that of BHM09al Thm 5.4], which is in turn 
largely inspired by that of IINor86l Thm 3.2]. We therefore only give a sketch of the 
proof and refer to IINor86l |HM09al INua06ll for details. Let ^ G R"* be an arbitrary 
deterministic vector with ||^|| = 1 and fix a terminal time t > 0. We then define the 
R^''^ -valued stochastic process by 

Note that t is 11 fixed terminal time, so that /i^ is a stochastic process which is not 
adapted to the filtration generated by the Wiener processes driving ( 14.1b . With this 
notation at hand, it follows from the integration by parts formula on Wiener space 
IINua06l that one has the identity 

t 



{DVMx),0='E(^^ixt)J h^{s)dW{s) 



where the stochastic integral should be interpreted as a Skorokhod integral IINua06ll . It 
then follows from the generalised Ito isometry IINua06ll that 

{DVMx),S,Y- <{Vtip^){x)^(^j^ \\h^(s)fds + \\%h^{s)fdsdrY 

where Si^F denotes the Malliavin derivative at time r of the random variable F. It 
follows from the chain rule that one has the identity 

Jo 

Note now that the Malliavin derivative of the Jacobian J^ ^ is given by 

{0 ifr>t, 
jfl^Ts, Js^rO Otherwise. 

Collecting all of these identities, it then follows immediately from the a priori bounds 
of Proposition |43] that the requested bound holds. □ 
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With all of these preliminary results at hand, we can now proceed to the proof of 
Theorem l4.4l 

Proof of Theorem \4.4\ As a first step, we show that the assumptions of Theorem 13.21 
hold. It follows from Proposition 14.71 that Assumption 4 is satisfied with Ui{x) = 
and U2{x) = exp((5||a;|p) for arbitrary choice of 5 > {). In view of Proposition |43] it 
is therefore natural to choose V{x) = W{x) ~ exp(77||x|p) for some sufficiently small 
value of 77 > 0. It is possible to first choose -q small enough so that the bounds in 
Proposition |43] hold and then 6 small enough so that the bounds ( 13.21 ) hold. 

We now want to show that ( I3.3l l holds as well. It follows from Assumption |7] 
combined with the support theorem IISV79II that Vt{x, B^{y)) > for every x,y £ R'', 
t > and £ > 0. Since furthermore the Markov semigroup Vt is Feller and the level 
sets of V are compact, ( 13.31 ) follows at once. We thus conclude that, for every t > 0, 
the assumptions of Theorem 13.21 hold, so that the semigroup Vt possesses a spectral 
gap in the space Cy y for sufficiently small values 77 > 0. 

We now show that Theorem l2.3l can be applied to our situation, if we choose a to 
be any of the components of either M or E. First of all, it follows from our previous 
discussion that Assumption [1] is satisfied for every t > 0, provided that we choose 
V{x) = 1 + /3exp(77||.T|[^) and W(x) = exp(?7||a;p) for sufficiently small values 
/?, 5, and 77 (possibly depending on the value of t). Furthermore, if we set U {x) = 
2exp(77'||a;p) for some sufficiently small value 7/ > 77, then it follows from (14.4b that 
Assumption|3]is also satisfied. It therefore suffices to check that Assumption |2] is also 
satisfied. 

Denote by $t(a;, M, S) the solution to ( 14.1b with initial condition x and parameters 
Af and E. By considering and E as additional dynamical variables that follow the 
evolution equations TV" = 0, E = 0, it follows from iKun90l ICE83l that on R'' x J^, 
$t is a C°° map of all of its arguments. Furthermore, it follows from the variation of 
constants formula that its derivative with respect to M and E in arbitrary directions SAf 
and (5E is given by 



Combining this with the bounds from Proposition l4.5l it follows at once that Assump- 
tion |2]is satisfied with 



provided that /3, 6, and 77 are sufficiently small. This allows us to apply Theorem l2.3l 
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